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Abstract 

A formulation for stationary axisymmetric electromagnetic fields in 
general relativity is derived by casting them into the form of an anisotropic 
fluid. Several simplifications of the formalism are carried out in order to 
analyze different features of the fields, such as the derivation of electromag- 
netic sources for the Maxwell field in the form of thin layers, construction 
of new solutions, and generation techniques. 
PACS: 04.20.Cv, 04.20. Jb, 04.40. +c 

1 Introduction 

The study of coupled stationary axisymmetric electromagnetic and gravitational 
fields is of great interest in astrophysics in order to provide a description for the 
exterior of steadily rotating objects. Although the whole system of equations 
can be reduced to two complex partial differential equations [1] and special 
techniques have been developed to generate solutions, for instance the ones 
reviewed in [2], [3] for vacuum fields, we still lack exact solutions whose physical 
interpretation is appealing. 

In section 2 we shall extend the exterior differential system formalism pre- 
sented in [4] and [5] for perfect fluids to electromagnetic fields. The Maxwell 
field will be presented as an anisotropic fluid. Null fields will not be considered 
since they have a null Killing vector [6] and have already been studied. The 
main advantage of using differential forms is that we are not forced to choose 
coordinates from the beginning and therefore there is more freedom for solving 
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the equations and simplifying the system in different ways. In section 3 a sim- 
plification will be carried out on this Maxwell fluid to cast it into a shear-free 
form. Advantage will be taken from it to construct magnetic dipole thin layers 
as surface sources for the magnetic field of static (among other) electrovacs and 
to extend previous results on angular momentum densities for vacuum metrics 
to some electrovacs. Section 4 is devoted to another simplification of the formu- 
lation. This time the Maxwell fluid will take the form of an irrotational fluid. 
An extension of some results obtained in canonical coordinates is provided. In 
section 5 it is performed a simplification of different nature: The Ricci tensor is 
diagonalized and a family of new electrovacs is constructed. 

2 Exterior differential system for stationary ax- 
isymmetric electrovacs 

2.1 Cartan's first structure equations 

We shall follow the approach developed in [4] and [5] for the kinematical quan- 
tities and define an orthonormal vierbein {9°, 9 1 , 9 2 , 9 3 } where {9°, 9 1 } lie on 
the space spanned by the orbits of the isometries {d t , whereas {9 2 , 9 3 } will 
lie on the orthogonal space. In order to partially diagonalize the electromagnetic 
stress tensor we shall choose 9° parallel to the direction of its timelike eigenvec- 
tor as the timelike leg of the tetrad. The four eigenvalues of the stress tensor 
have the same absolute value and one of them is negative. This one corresponds 
generically to the eigenvector which lies in the direction of the magnetic and 
electric field in a frame in which they are parallel or, if the Lorentz invariant 
E ■ B is zero, then the respective eigenvector lies in the direction of the elec- 
tric (magnetic) field in a frame in which the magnetic (electric) field vanishes. 
Therefore we shall be able to view the electromagnetic field as an anisotropic 
fluid with the absolute value of the pressure equal to the density and consider 
u = —9° as its velocity form. Imposing that the tetrad be torsion- free we get 
the following Cartan equations: 



where a is the acceleration of the 'fluid', w — *lo {ui is the vorticity form derived 
from u and * stands for the two-dimensional Hodge dual in the space orthogonal 
to the orbits of the Killing vectors {d t , c^}), s is a one-form associated with the 
shear of the fluid, b = d\np (p is the radial pscudocylindrical Weyl coordinate) 
and v is just a connection in the two-space spanned by 9 2 and 9 3 . We still have 
a great amount of freedom for further simplifications. 



du = aAu + wA9 1 

d9 1 = (b - a) A 9 1 + s A u 
d9 2 = —v A 9 3 
d9 3 = v A 9 2 



(1) 
(2) 
(3) 
(4) 
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2.2 Bianchi identities 

Bianchi identities are easily obtained by exterior differentiation of the previous 
set of equations: 

db = (5) 

da = w A s (6) 

dw = -(6 -2a) Aw (7) 

ds = (b - 2a) A s (8) 

2.3 Maxwell's equations 

Now we shall describe the Faraday and Maxwell electromagnetic strength two- 
forms (F and its four-dimensional dual 4 * F ) in terms of the electric and 
magnetic one-forms, respectively E and B) . In what follows we shall restrict 
ourselves to electromagnetic forms E and B with no projection on the orbits of 
the Killing vectors: 

F = -E Au + *B A6 1 (9) 
4 *F = B Au + *E A6 1 (10) 
Therefore, the Maxwell vacuum equations dF = and d 4 * F = imply: 

d£-£Aa + *5As = (11) 

d*5 + £Aw-*5A(6-a) = (12) 

d*E-BAw-*EA(b-a) = (13) 

dB-BAa-*EAs = (14) 

These equations are not a consequence of the integrability conditions for the 
Einstein equations unless the Lorentz scalar E ■ B is different from zero. These 
equations can be written in a more compact expression defining a complex one- 
form /: 

f = E + iB (15) 
df = —a A f — is A */ (16) 
d * / = (a - b) A */ + iw A f (17) 
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2.4 Einstein's equations 

Considering that the electromagnetic stress tensor has the following expression: 

Ta b = -^{ paCF bc - \5 a b F cd F cd }, (18) 
the Einstein equations take the general form: 

d * (w - s) + 2a A *w + 2 (a - b) A *s = AE A B = 2if A / (19) 

d*a + bA*a+ A *w — A *s = E A *E + B A *B = f A *f (20) 

d*b + bA*b = (21) 

db + b A b - ^ (s - w) A (s - w) + 2a A a - 2b A a + 

+2v A*b = 2E AE + 2B AB = 2f Af (22) 

d * b + b A *b — i(s — w) A *(§ — w) + 2a A *a — 2b A *d — 

—2v Ab = 2EA*E + 2BA*B = 2fA *f (23) 

dv + a A *b - a A *a + i(s - w) A *(s — w) = (24) 

The symbol ~ denotes the linear transformation in the {9 2 , 9 3 } -space defined 
by the following expression: 

§ 2 = 6 2 § 3 = -9 3 (25) 

As in the perfect fluid case, equations 22 to 24 can be left for the end, since 
the connection v can be algebraically obtained from them and then equation 24 
is automatically satisfied. 

3 The electromagnetic field as a rigidly rotating 
fluid 

3.1 Ernst's equations 

As we have already said, our exterior system can be simplified in many ways. 
This freedom comes from the possibility of describing the electromagnetic field in 
different reference frames preserving the structure of the equations. We can shift 
from one frame to another performing a local boost or rotation or a combination 
of both. 
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One boost that can always be applied is the one which takes us from our 
original frame to another in which the velocity u is shear- free (s = 0). This 
transformation is described by the following equations: 

u' = cosh A ill sinh A 6 1 (26) 

9 1 ' = sinhA w + coshA 9 1 (27) 

dX = cosh A sinh A (b' - 2a') + cosh 2 As' - sinh 2 A w 1 (28) 
The electromagnetic field changes in the following way: 

E = cosh A E' - sinh A * B' (29) 

B = cosh A B' + sinh A * E 1 (30) 

This transformation is compatible with the whole set of equations as it hap- 
pened in the vacuum case [4] . Therefore there is no loss of generality in taking 
,s = from the beginning. 

We can formally integrate equations 5 to 7 and 16 to get: 

a = dU b = d\np (31) 
w = p- x e 2V dA (32) 

/ = -e- u d$ (33) 

From equation 19 we obtain the exact differential that defines the twist 
potential x 

d x = e 2U * w + 2i$d$ (34) 
and construct the Ernst potential e in the following way: 

de = de 2U + id X Tie = e 2U - l>$ (35) 

Substituting 33 in 17 we get an elliptic equation for the scalar complex 
electromagnetic potential $: 

d*d$ + b/\*d<i> = — ^{de + 2i>rf$} A *d$ (36) 

Tie + 

And if we combine the previous equation with the Raychaudhuri equation 
(20), and 34 with 31 and 7, we obtain another second order elliptic equation, 
namely the Ernst equation: 

d*de + bh*de= — ={de + 2$d&} A *de (37) 
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The whole set of equations for stationary axisymmetric electrovacs reduces 
to these two complex equations for the potentials e and $ [1]. 

The equations 5 and 21 just introduce two functions p and z, the Weyl 
pseudocylindrical coordinates, that are related by the Hodge dual of their dif- 
ferentials: 

*dp = — dz (38) 

And therefore the metric takes the following form, after integrating the Car- 
tan equations: 

ds 2 = -e 2U (dt - Ad4>f + e- 2U [e 2k (dp 2 + dz 2 ) + p 2 d4> 2 ] (39) 

3.2 Construction of electromagnetic dipole surface densi- 
ties 

This formalism can be used to obtain electric and magnetic dipole surface den- 
sities for the sources of some asymptotically flat stationary axially symmetric 
solutions of the Einstein-Maxwell equations in a similar way as it was done in 
[7] to construct angular momentum densities for vacuum metrics. The formulae 
derived by Israel in [8] cannot be applied since they only yield non- vanishing ex- 
pressions for the mass, charge and angular momentum, but not for the magnetic 
momentum. 

We shall only consider static fields and also nonstatic metrics for which 
w A / — -if we want to calculate only magnetic (electric) moment densities, we 
shall only need w A E = (w A B = 0)-. By asymptotically flat we mean that 
the metric has the following behavior at infinity (r — > oo) in some coordinates 
(t,r,0,4>): 



, 9 /, 2m. , , 2Jsin 2 9 , , x9 

ds 2 = -(1 ){dt + d<j)) 2 + 

r r 

+ (1 + — )[dr 2 + r 2 (d0 2 + sin 2 0# 2 )] (40) 

r 

where m is the source's total mass and J is the total angular momentum. There- 
fore the metric functions e 2U and A and the twist potential \ (since A and \ 
are related by the dual of the vorticity one-form) take this form at infinity: 

^!-^ + 0(l) A.-^l^) («) 
X-^+Oa) (42, 
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To calculate the magnetic dipole we shall require that the magnetic scalar 
potential, that is the imaginary part V of the potential $ must be asymptotically 
dipolar: 

V^ + 0(1) ,43, 

As it happens in flat spacetime, the magnetic scalar potential can be defined 
only outside the sources. If we have a sheet of dipoles as a source, the limit 
values of the potential on either side of the surface will be different and this 
discontinuity reveals the existence of the source (cfr. for instance [9]). Our 
purpose will be to extend this result to curved spacetimes. 

Following the approach developed in [7], we integrate equation 12 with the 
previously stated restrictions, that combined with 33 provides us with two ex- 
pressions for the magnetic field: 

B = X{f} = -e~ u dV = -p- 1 e u *dW (44) 

V = (45) 

This function W has the following asymptotic expression due to the behavior 
of V at infinity: 

Msin 2 6» /1\ 
W=^~ + [y,) ( 46 ) 

And now we shall integrate the scalar product of the magnetic field B and 
the differential of a function Z(r, 9) to be determined. The domain of integration 
will be the space V 3 orthogonal to the congruence defined by u. The metric on 
this space is the projection g = 4 g + u®u. Using both expressions for B we get: 

= / y/g < [B + *(*B)],dZ > dx 1 dx 2 dx z = 
Jv 3 

= I ^{-e- u g^d„Vd v Z + p~ 1 e u e^d„Wd u Z}dx 1 dx 2 dx 3 (47) 
Jv 3 

In order to express the integrand as a total derivative we choose Z so that 
it fulfills the following differential equation: 

d„{Jge- u g» v d v Z) = G (48) 

As an asymptotic boundary condition we impose that Z behaves like r cos 9 
at infinity. 

The potential V is discontinuous across a closed surface S if we are to have a 
magnetic moment surface density, therefore we have to split the space into two 
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pieces V 3 + and V^ - , respectively the outer and inner part of V3. If the surface 
S is open, it can be extended to a closed one, taking the value zero for the 
discontinuity on the closure. V + and V~ shall denote the limit values of V on 
either side of the surface. 

Since the Levi-Civita tensor on the space orthogonal to the orbits of the 
Killing fields is = e 2 ( u ~ k ) [^] the integrand in 47 can be written as: 

d»{-^e- u g^Vd v Z + W[^]d v Z} (49) 

We have the desired total derivative and we can express the integral as 
a surface integral on the boundary £ = dV$~ U dV^ with unitary normal n. 
The boundary of consists of S and the sphere at infinity S 2 (oo) and the 
boundary of is just S. Taking into account the required behavior of the 
metric functions and the magnetic potential at infinity, the integral on S 2 (oo) 
yields: 



/ dS[V]e~ u g^n„d„Z = 4ttM (50) 
Js 

And therefore we obtain the following formula relating the jump of the mag- 
netic potential [V] across S with the total magnetic moment M: 

AttM = ( dS g (51) 
Js 

a= ^[V]e- u g^n^Z (52) 

47T 

We can consequently infer that a is the source's magnetic moment density 
for the Maxwell field. 

Similar calculations can be done with the real part of the electromagnetic 
potential <J> to compute electric dipole surface densities. 



3.3 Angular momentum surface densities 

The formalism developed in [7] for interpreting discontinuities of the twist po- 
tential of asymptotically flat stationary axially symmetric vacuum metrics can 
be extended in a straightforward way to cope with electrovac solutions with 
E A B = 0, since the only equations involved are the Bianchi equations and the 
equation for the vorticity, 19, and these remain the same as in the vacuum case 
if we impose the mentioned restriction. Hence the angular momentum surface 
density for stationary axisymmetric vacuum metrics and electrovacs with paral- 
lel magnetic and electric fields has the following expression in terms of the jump 
[x] of the twist potential: 

<Jrot^—^[x]e- ZU 9^n^d u Z rot (53) 
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where the function Z ro t is asymptotically r cos 9 and satisfies the following dif- 
ferential equation: 

d^e- 3U g^d v Z rot ) = (54) 
3.4 Bonnor's transformation 

One of the peculiarities of this formalism for constructing surface densities is its 
behavior under Bonnor transformations. These transformations are a method 
for generating magnetostatic electrovacs from stationary nonstatic metrics [10]. 
As the solutions of the Einstein-Maxwell equations arc fully characterized by 
the Ernst and electromagnetic potentials, it will suffice to give the rules that 
determine the potentials for the electrovac: 

£ &vac &vac (^5) 

V = ixvac (56) 

In order to have a real solution the magnetic potential has to be real, there- 
fore the parameters have to be rearranged with a complex transformation to 
achieve this purpose. 

From these expressions it is clear that e u = e^ c . Since Z is a function only 
of p and z, in canonical coordinates, the only components of the three- metric g 
of V3 that appear in the differential equation 48 are g pp — g zz — e 2U ~ 2k . The 
square root of the determinant of this metric is ^fg = e~ 3U+2k p and therefore 
equation 48 takes this form in terms of the original vacuum metric functions 
after rearranging the parameters: 



d^y/ge-VgTdvZ) = c^e"^ pd^Z) = 

d^(e-^p5^d v Z) = d^^- c (e- 3U ) vac gZ c d u Z) = (57) 

That is, the function Z needed for the construction of the magnetic dipolc 
surface density satisfies the same differential equation as the Z rot involved in 
the calculation of the angular momentum density for the original vacuum metric 
g vac . Therefore the same function is valid for both solutions, once we have 
reinterpreted the parameters, and the formalism is compatible with the Bonnor 
transformation. 

3.5 An example: Bonnor's massive magnetic dipole 

As an example of how this formalism works, we shall calculate the magnetic 
source for Bonnor's massive magnetic dipole [11]. This is a Bonnor transforma- 
tion of the Kerr metric and therefore we can use the Z function obtained in [7] 
for its source's angular momentum density. 
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Bonnor's metric has the following form: 



ds* = -(i- 2 2 7 2 j 2 dt a + 

r z — cos^ # 

+ (1 - r2 _ 2 Q T c r og2 , )- 2 {(r 2 - a 2 - 2™r) sin 2 ^ 2 + 
(r 2 - a 2 cos 2 0-2™^ (rfg2 + rfr 2 j 



[(r — m) 2 — (a 2 + to 2 ) cos 2 #] 3 r 2 — 2mr — a 

from which we can read the necessary metric functions: 



\ 2 

2mr \ 4to „ / 1 



e = 2-^TT =1 + - (59) 

\ r z — a A cos A 9 / r \r z ) 

The magnetic field is determined by the scalar potential: 

2amcos9 2amcos8 n ( 1 \ 

V = - 2 = Ya = 5 + -4 60 

r A — cr cos A 9 r z \ r J 

Therefore both the gravitational and magnetic fields fulfill the required 
asymptotic conditions with mass equal to 2m and magnetic moment 2am. 

If we consider the metric to be written in a sort of oblate spheroidal coordi- 
nates as it is done in [8] for the Kerr metric, then events on the surface r = 
with polar angle 9 have to be identified with events with polar angle equal to 
it — 9. To avoid double-counting these points, the 9 coordinate will range from 
to tt/2 on this surface. Therefore the function cos# suffers a discontinuity upon 
crossing r = 0, since it shifts sign from positive in the upper subspace (r > 0, 
< 9 < n/2) to negative in the lower subspace ( r > 0, tt/2 < 9 < tt). 

Hence the magnetic potential V = !{$} is discontinuous on the surface 
r = 0. The difference between the values taken on the upper and the lower side 
of the surface is: 

a cos 9 

The Z function satisfying 48 can be obtained from the Z rot for Kerr, as it 
was shown in the previous section, taking into account that Hence, 

„ 2a 2 TO(5cos 3 6» - 3cos6>) , N 

Z=(r- 3m) cos 9 - — 62 

5(r J — a 1 cos 2 - 9) 

The surface element for r = is: 

5 

\d cos 9 — to sin 9\ t3 ' z 

We are ready now to write the magnetic moment surface density for the 
source of Bonnor's magnetic dipole using the formula 52: 
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m | a 2 cos 2 9 - m 2 sin 2 6 | 3/2 

^4 COS 4 # 



(64) 



The integral defining the total magnetic moment yields the expected result: 

(■7T/2 



/■ Z7T Pit I £ 

= I a (IS = I / d9 (ma sin 0) = 2ma (65) 
Jo Jo 

But unfortunately the source lies in a region where the signature of the 
metric is not the usual one (for instance, the angle (p is no longer a spacelike 
coordinate). Any other surface source would have to include the region S in 
order to cope with the discontinuities of the magnetic potential. Therefore 
its interpretation as a physical material source is dubious and one would have 
to resort to three-dimensional sources to hide the unphysical regions of the 
clcctrovac metric. 



4 The electromagnetic field as an irrotational 
fluid 

4.1 Ernst's equations 

Instead of performing a boost in order to cancel the shear one-form s, as it was 
done in the previous section, we can always write our equations in a frame where 
the dual of the vorticity w is zero. This is also compatible with our system of 
equations, as it happened for vacuum fields [4], [5]. Of course, this would not 
be possible in general for a perfect fluid. The differential equation defining the 
A parameter of the boost is now: 

dX = cosh A sinh A (2a' - b') + cosh 2 A w' - sinh 2 A s' (66) 

This amounts to many simplifications, as it happened in the shear-free for- 
malism. The equations 5 to 8 and 17 can be integrated to get: 

a = dU b = d\np (67) 
s = pe- 2U dC (68) 

*/ = -e u p~ 1 d^ (69) 

From equation 19 we get a exact differential that can be used to define what 
could be called 'deformation potential' ip in analogy with the twist potential: 

d?/> = p 2 e- 2U * s - 2i*d* (70) 
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and construct another complex potential r\ similar to the Ernst potential: 



dr) = d{-p 2 e~ 2U ) - idip 1Z{ri} = -p 2 e~ 2U - (71) 

If we substitute 69 in 16 another generalization of the Laplace equation is 
achieved, this time for the potential 

d * + b A *d^ = -{dn + 2*d*} A *d* (72) 

TZr] + $$ 

The analog of the Ernst equation can be obtained in this irrotational formal- 
ism in a similar way taking into account the latter equation, the Raychaudhuri 
equation (20) and the equations 70, 67 and 8: 

d* d-q + b A*di] = — ^ ^^ {drj + 2^d^} A *drj (73) 

These two equations, together with the coordinate condition *dp = dz, com- 
prise the whole problem of calculating stationary axially symmetric electrovac 
solutions. In this formalism the metric takes a form different from the canonical 
one: 

ds 2 = ~e 2U dt 2 + e- 2U [e 2k (dp 2 + dz 2 ) + p 2 (d<P + Cdt) 2 ] (74) 

Both sets of equations, the Ernst equations (36 and 37) and these new Ernst- 
like equations, look the same exchanging e for r\ and $ for \1/ as it happened in 
the vacuum case [4], [5]. However, there is a difference that did not occur for 
vacui and that is hidden behind the minus sign that it is written before the norm 
of the rotation Killing vector, ||c^|| = p 2 e~ 2U , as it appears in the definition for 
the 77 potential. This sign enforces the real part of rj to be negative if it has to 
deal with a physical situation. Otherwise the rotation Killing vector would be 
timelike. On the other hand e can be either positive or negative. Besides, the 
conditions on the potentials for the metric to be asymptotically flat are rather 
different: 

r, = -r 2 sin 2 6 + 0(r) (76) 

Although the equations have been cast in the same form, in fact the poten- 
tials have very different meaning. 

One could think of the following transformation in order to generate new 
stationary axially symmetric electrovacs as the electromagnetic generalization 
of the transformation introduced in [4] and [5]: 

e — ► V (77) 
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$ — > * (78) 

But this would lead to unphysical metrics unless the seed metric has (H ) 

signature, since this transformation induces the following exchange between the 
norms of the Killing vectors: 

e 2U — > -p 2 e- 2U (79) 

It could be used then to generate new solutions from the forbidden regions 
of other spacetimes. 

4.2 Two families of 'irrotational' solutions 

The first attempt for constructing new solutions of the Einstein-Maxwell equa- 
tions with the irrotational formalism could be to make use of the ansatze that 
have been profitable for solving the standard Ernst equation. For a start it can 
be shown that the Weyl family of static electrovacs [12] can be extended to yield 
new solutions. 

The Weyl ansatz, in terms of the Ernst potentials, assumes a functional 
dependence between e and $, that now are real (Since the metric is static, the 
twist potential can be taken to be null): 

e = e($) (80) 

This ansatz leads to two equations, that can be integrated, considering the 
asymptotic flatness condition for the Ernst potential (75): 

e(*) = e 2U = 1 - 2c$ + $ 2 (81) 

where c is a constant. 

d * d$ + b A *d$ = ~ 2c + 2$ d<5> A (82) 
1 - 2c$ + $ 2 v ' 

that can be solved in terms of a function Y, solution of the reduced Laplace 
equation: 

d * dY + b A *dY = (83) 

via the following relation: 

f -V^^cothF + c f (c 2 - 1) sinlT 2 Y c 2 > 1 

$ = \ - y 1 +c e 2U = 1 Y 2 c 2 = 1 (84) 

{ -VI - ^coty + c { (i-c 2 )shy 2 y c 2 < 1 

A similar ansatz can be applied in the irrotational formulation of static 
axisymmetric electrovacs: 
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V = V&) (85) 

The equations that this ansatz yields are entirely similar to those of the Weyl 
ansatz: 

) = p 2 e- 2U =c' + 2c* - * 2 (86) 

d * d* + 6 A *<i* = , 2 ^ TO <i* A *<i* (87) 

c' + 2c* - * 2 y ' 

For physical reasons we want c' + 2c* — * 2 to be positive and this can only 
be achieved if d + c 2 > 0. With this restriction in mind the physical cases to 
be considered reduce to just one: 

* = c + yjd + c 2 tanh Y p 2 e~ 2U = (c' + c 2 ) cosh" 2 Y (88) 

where Y is again a solution of the reduced Laplace equation. 

From the previous equation we learn that p 2 e~ 2U is uppcrly bounded by 
d + c 2 and hence these solutions cannot be asymptotically flat at infinity, since 
the norm of the rotational Killing vector is expected to diverge at infinity. This 
family could be used to represent an electrovac spacetime surrounded by a ma- 
terial source to which it should be matched. 

Another ansatz that can be extended to this irrotational formalism is that 
of Papapctrou [13] for vacuum metrics. Following [14] we can write this ansatz 
as an assumption of functional dependence between the real and the imaginary 
parts of the Ernst potential: 

e = e 2U +i X U = U( X ) (89) 
The real and imaginary parts of the Ernst equation imply: 



e 



iU = c' + 2c X -x 2 (90) 



where c and d are constants. 



d*d X + b/\*d X = , 2C 2X 2 d X f\*d X (91) 
d + 2c X - x 2 

which is similar to the one we have integrated for the extension of the Weyl 
solution (87). Therefore the twist potential has the following expression in 
terms of an arbitrary solution of the reduced Laplace equation: 

X = c+ y/d + c 2 tanhY" (92) 

The restriction c' + c 2 > has to be imposed again for physical reasons. The 
Ernst potential takes the form: 
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e = y/d + c 2 (cosh l Y + i tanhF) + ic (93) 

If we apply this ansatz to the irrotational formalism, p 2 e~ 2U = g(ip), instead 
of to the habitual shear-free formalism, then the equations have the same form 
and therefore we have a different family of solutions: 

T] = \fc!~~+~(? (coshT 1 Y + i tanh Y) + ic (94) 

The same considerations for the asymptotic behavior of the extension of the 
Weyl family apply to this extension of the Papapetrou class. 

5 Diagonalization of the stress tensor 

Classical nonnull electromagnetics fields can be put in a frame where their stress 
tensor is diagonal by performing a Lorentz transformation. It implies taking 
the Poynting vector to be zero by letting the electric and magnetic field become 
parallel (if the Lorentz invariant E ■ B is different from zero) or by cancelling 
either the electric or the magnetic field (if E ■ B = and the other Lorentz 
invariant ||-E|| 2 — ||-B|| 2 is respectively negative or positive). 

The same simplification can be applied without loss of generality in general 
relativity. This amounts to taking EAB = in the exterior system and therefore 
the Einstein equations for the difference of the vorticity and the shear and for 
the dual of the acceleration remain as follows: 

d* (w - s) + 2a A *w + 2 (a - b) A *s = (95) 
11 

d*a + bA*a+ -w A *w — -s A *s = E A *E + B A *B = f A */ (96) 
5.1 'Seminull' Maxwell fields 

In this section we shall call 'seminull' the Maxwell fields that have one of the 
Lorentz invariants identically zero in opposition to null and generic fields. There 
are three different cases: 

As we have already mentioned, when E ■ B = and E 2 — B 2 < the 
equations can be expressed in a frame where E = 0. We shall write only the 
relevant equations where the electromagnetic field appears: The Raychaudhuri 
equation (20) and the Maxwell equations: 

d*a + bA*a+ A *w - A *s = B A *B (97) 
dB - B A a = (98) 
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d* B - *B A (b- a) = 



(99) 



B A w = (100) 

*BAs = (101) 

From the last two equations we get as a constraint that the vorticity and 
shear one-forms have to be orthogonal. 

A similar case happens when E ■ B — and E 2 — B 2 > 0. Now it is B 
the field that can be taken to be zero. The resulting equations are the same as 
97-101 after changing B for E. 

The last case of seminull fields consists of those that have E ■ B ^ 0, but 
E 2 = B 2 . Hence, in the frame where both fields are parallel, we have E = ±B. 
Therefore E and B fulfill the same Maxwell equations, which take the form of 
equations 98 to 101. The Raychaudhuri equation becomes slightly different: 

d*a + b/\*a+ A *w — A *s — 2E A *E (102) 

From the whole set of equations for the seminull Maxwell fields we get as a 
consequence that a solution for one of the cases provides also a solution for the 
other two by the following transformation rule: 

Be?<b 2 — * E E 2 >B 2 — ► -E E -B=n = ±2^ EB=0 (103) 

Although the three cases are physically different, their mathematical de- 
scription is the same. 

5.2 A family of seminull electrovacs 

In order to derive new solutions of the Einstein-Maxwell equations we shall make 
use of an ansatz for the 'magnetic seminull' field: 

We take it to be shear-free (s — 0) and assume the following condition on 
the kincmatical one-forms: 

^wA*w = BA*5 (104) 

Due to the constraint 100, this ansatz implies that B = ±^w. Therefore 
we have two different sets of equations for the vorticity and its dual: 

dw = -(6 -2a) Aw (105) 
d*w + 2aA*w = (106) 
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dw — w A a = 



(107) 



d* w - *w A (b - a) = (108) 
that will introduce a new constraint: 

6 = 3a (109) 

The integrability conditions for this constraint impose no new restriction 
since now the equations for the acceleration form and its dual are consistent 
with those governing b. 

We already know that the equations for the one-form b, 5 and 21, just in- 
troduce the Weyl coordinates p and z, such that *dp — dz. Hence, the only 
equations that remain to be solved are 105 and 106, that can be straightfor- 
wardly integrated: 

w = p~idA (110) 

*w = p~^dx (HI) 
whose integrability condition in Weyl coordinates is: 

A pp + ^-A p + A zz =0 (112) 
op 

which is a well-known equation in hydrodynamics, the Tricomi equation in 
canonical coordinates when the parameters allow to classify it as an elliptic 
equation. 

The twist potential can be obtained from the solutions of the Tricomi equa- 
tion by quadratures: 

Xp = P*A x (113) 

Xz = -P*A P (114) 

Since a = |c?lnp, the metric can be expressed as follows in canonical coor- 
dinates: 

ds 2 = -cpi (dt - Ad<P) + c^p-i [p 2 d<p 2 + e 2k (dp 2 + dz 2 )} (115) 

c being a positive constant and k is obtained from the Cartan equations once 
the connection v is calculated from the equations 22 and 23. If we write the 
tetrad forms Q 2 and 9 3 as: 

9 2 =e Q dz 8 3 = e Qdp (116) 
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then the metric function Q = k — U can be integrated from this set of equations: 

Qz = -\p^A p A z (117) 

Q P = -|p- 1 -|p*(V-V) (H8) 

whose integrability condition is precisely the Tricomi equation 112. 

From the expression for the dual of the magnetic one-form *B we get the 
Faraday two-form F: 

F = ±-^=dx A d<p (119) 
V2c 

and the magnetic field is : 

B = ±-^e- Q p-i(A p 9 2 +A Z 9 3 ) (120) 
V2 

From the general expression for the metric it is obvious that it is not asymp- 
totically flat in canonical coordinates. This seems to arise from the fact that the 
gravitational field obtained from the goo component has cylindrical symmetry 
and therefore cannot be due to a compact source. The magnetic field, however, 
is not cylindrical unless we assume that the metric function A does not depend 
on the coordinate z, in which case we would have an extra Killing vector d z - 

This family of solutions of the Einstein-Maxwell equations has got generically 
only two Killing vectors and Petrov type I. It is also generically nonstatic. 



6 Discussion 

In this paper it has been introduced a new formalism for studying stationary 
axially symmetric electrovacs making use of exterior differential systems of equa- 
tions and putting the electromagnetic field into the form of an anisotropic fluid. 
This has been useful since we have not attached ourselves to any special coor- 
dinate system from the beginning and therefore we have been able to simplify 
our equations in different manners according to our needs. 

This formulation of the coupled gravitational and electromagnetic fields has 
given results both in the interpretation and in the derivation of solutions of 
the Einstein-Maxwell equations. It has served to deal with sheets of dipoles as 
sources for the fields and, on the other hand, to extend old families of electrovacs 
and also to construct new ones. In particular, the result concerning the magnetic 
dipole sources for the fields is of great interest since they cannot be obtained by 
the usual formulation derived by Israel [8]. 

For a future work it would be interesting to extend further the applicability 
of the techniques developed so far for interpreting the discontinuities of the 
Ernst potentials so that other physical situations fit in the formalism. This 
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would provide a deeper insight into general relativity and its comparison with 
flat-spacetime Physics. 

The present work has been supported in part by DGICYT Project PB89-014-2; 
L.F.J, is supported by a FPI Predoctoral Scholarship from Ministerio de Edu- 
cacion y Ciencia (Spain). The authors wish to thank L.M. Gonzalez- Romero 
and J. A. Ruiz-Martin for valuable discussions. 

References 

[1] F. J. Ernst, Phys.Rev. 168, 1415 (1968) 

[2] Solutions of Einstein's equations: Techniques and Results (eds.: C. Hoense- 
lacrs and W. Dietz, Springer Verlag), Berlin-New York (1984) 

[3] El Escorial Summer School on Gravitation and General Relativity 1992: 
Rotating Objects and Other Topics (eds.: F. J. Chinea and L. M. Gonzalez- 
Romero), Springer- Verlag, Berlin-New York (1993) 

[4] F. J. Chinea and L.M. Gonzalez-Romero, Class. Quantum Grav. 9, 1271 
(1992) 

[5] L. M. Gonzalez-Romero, Ph.D. thesis, Universidad Complutense de Madrid 
(1991) (unpublished) 

[6] M. Gurses, J. Math. Phys. 18, 2356 (1977) 

[7] L. Fernandez- Jambrina and F.J. Chinea, Phys. Rev. Lett. 71 2521 (1993) 

[8] W. Israel, Phys. Rev. D 2, 641 (1970) 

[9] O. D. Kellogg, Foundations of Potential Theory. Dover, Nueva York. 

[10] W. B. Bonnor, Z. Phys. 161, 439 (1961) 

[11] W. B. Bonnor, Z. Phys. 190, 444 (1966) 

[12] H. Weyl, Ann. d. Phys. 54, 117 (1917) 

[13] A. Papapetrou, Ann. d. Phys. 12, 309 (1953) 

[14] D. Kramer, H. Stephani, M. MacCallum, E. Herlt, Exact Solution's of Ein- 
stein's Field equations Cambridge University Press, Cambridge (1980) 



19 



